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Abstract. We formulate the dual fermion approach for strongly correlated electronic 
systems in terms of the lattice and dual effective interactions, obtained by using 
the covariation splitting formula. This allows us to consider the effect of six-point 
one-particle reducible interactions, which are usually neglected by the dual fermion 
approach. We show that the consideration of one-particle reducible six-point (as well 
as higher order) vertices is crucially important for the diagrammatic consistency of 
this approach. In particular, the relation between the dual and lattice self-energy, 
derived in the dual fermion approach, implicitly accounts for the effect of the diagrams, 
c"j ' containing 6-point and higher order local one-particle reducible vertices, and should be 

applied with caution, if these vertices are neglected. Apart from that, the treatment 
of the self-energy feedback is also modified by 6-point and higher order vertices; these 
vertices are also important to account for some non-local corrections to the lattice 
self-energy, which have the same order in the local 4-point vertices, as the diagrams 
| usually considered in the approach. These observations enlighten an importance of 6- 

point and higher order vertices in the dual fermion approach, and call for development 
of new schemes of treatment of non-local fluctuations, which are based on one-particle 
irreducible quantities. 
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1. Introduction 

Strongly correlated electron systems became one of the touchstone of modern 
physics. They demonstrate a variety of phenomena: magnetism, (unconventional) 
superconductivity, "colossal" magnetoresistance, and quantum critical behavior. The 
dynamical mean-field theory (DMFT)pQ E] allowed to describe accurately the Mott- 
Hubbard metal-insulator transition. [3]. DMFT becomes exact in the limit of high spatial 
dimensions (d — )> oo) and accounts for an important local part of electronic correlations. 
Real physical systems are however one-, two-, or three-dimensional. Therefore, nonlocal 
correlations, which are neglected in DMFT, may be important. Recently, a progress 
to go beyond DMFT through cluster extensions [U [5j El [7J [8] was achieved. These 
correlations are however necessarily short-range in nature due to numerical limitations 
of the cluster size [9]. 

This limitation motivated developing the diagrammatic extensions of the dynamical 
mean-field theory. The dynamical vertex approximation (DrA) was introduced in Refs. 
[TU| [TT] [T2| [T3| [TH H5] . Starting from the local particle- hole irreducible vertex, this 
approximation sums ladder diagrams for the vertex in the particle-hole channel, where 
particle-hole irreducible vertices assumed to be local, but the effect of the non-locality 
of the Green functions is considered. Alternative dual fermion (DF) approach was 
proposed in Refs. p2J [TTJ UHl EE] , which splits the degrees of freedom into the local ones, 
treated within DMFT, and the non-local (dual) degrees, considered perturbatively, with 
a possibility of summation of infinite series of diagrams for dual fermions [HI [20] . 

Although both abovementioned approaches use 4-point local vertex as an effective 
interaction between fermionic degrees of freedom (lattice fermions in case of DTA and 
dual fermions in the DF approach), they make in fact very different assumptions on the 
neglect of higher-order local vertices. Indeed, DTA operates with one-particle irreducible 
(1PI) vertices, and neglects six-point and higher order 1PI local vertices. At the same 
time, DF representation does not use the assumption of the one-particle irreducibility; 
in particular its formulation in Refs. [HI [TFJ HH [19] neglects one-particle reducible 
six-point and higher vertices. 

This difference appears to be important for analysing diagrammatic consistency 
of the abovediscussed approaches. While the dynamic vertex approximation is based 
on the diagrammatic approach, formulated in terms of the original lattice degrees of 
freedom, the diagrammatic consistency of the dual fermion approach (in terms of the 
same lattice degrees of freedom) has to be verified. In the present paper we show that 
the inculision of the one-particle reducible six-point (and more generally, higher vertices) 
into the DF approach appears to be necessary to make it diagrammatically consistent. 
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2. The model, dynamical mean-field theory, and the dual fermion approach 

2.1. The model and dynamical mean- field theory 

We consider general one-band model of fermions, interacting via local interaction 

H = ^ £ KvK,aC*,a + ^ Hint [Cia , C+ ] , ( 1 ) 

k,cr i,cr 

where Ci a ,cl a are the fermionic operators, and Ck tU ,c^ a are the corresponding Fourier 
transformed objects, a =t, 4- corresponds to a spin index. The model is characterized 
by the generating functional 

Z[ri,r] + }= J d[c,c + }exp{-S[c,c + ]+r] + c + c + ri} (2) 



S[c,c+] = |rfr|^4(r)^c 4(T (r) + i7[c,c + ]| 



(3) 



where c ia , cf a , r] icT , r]f a are the Grassman fields, the fields r] i(T , rif a correspond to source 
terms. The dynamical mean-field theory corresponds to considering the effective single- 
site problem with the action 

<Wt[c, c + ] = yidr{ -C\r) V i(r)q ff (r) + V H ini [c ia , c+] J> , (4) 



i (. a a 



where the "Weiss field" function C(r) and its Fourier transform C(^n) has to be 
determined self-consistently from the condition 

GiocOn) = . u . v 1 v ,. v = V G(k, zV n ) (5) 
C 1 {iv n ) - S loc (zz/ n ) ^ 

where 

G(k,iu n ) = G k = [G~£ - S loc (iz/ n )] _1 (6) 

G ok = i u n — £k is the lattice noninteracting Green function (we use the 4- vector notation 
k = (k, iv n )) and Si oc (w n ) is the self-energy of the impurity problem (0J, which is in 
practice obtained within one of the impurity solvers: exact diagonalization, quantum 
Monte-Carlo, etc. 

2.2. The formulation of the dual fermion approach by means of covariation splitting 
formula 

The dual fermion approach of Refs. [161 EES US] can be conveniently formulated in 
terms of an effective interaction of the lattice theory (see, e.g. Ref. |21] ) 

V[V, = - In / d[c, c + ] exp I £ G ok { C L + vL) ( c k* + Vk*) 

J I k,a 

-H int [c,c + }} 

= - In Z\G~lr] ka , Goiter] - v£a G okVk*- (7) 
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Expansion of the effective interaction V[r/,?7 + ] in source fields generates connected 
(in general, one-particle reducible) Green functions, amputated by the non-interacting 
Green functions of the lattice theory Gok- The relation between one-particle reducible 
and 1PI counterparts of the vertices can be involved. In particular, the (one-particle 
irreducible) self-energy of the lattice problem (i.e. the 1PI 2-point vertex function) 
can be extracted from the two-point connected vertex function via the relation 

v£ 2 = (i — CofcSfe). 

To split the local and non-local degrees of freedom in the effective interaction (CO) 
we use the covariation splitting formula [21], which is based on the identity 



(8) 



Yl C k lc ka°ka = In J d[c, c + ] exp I c t* B k 

k,cr " (. k.cr 

ka ) 

with Ak + Bk = Ck, Eq. (JSj) can be proven by integrating over the c,c + fields. For 
A k = ((iv n ), C k = G ok this implies 

V[t7,77 + ] = -In J d[c,c + ] exp l^cl a G- k % <a -V DM ft[v + + c + ,V + ^j (9) 

where Vdmft[?7, V + ] is an effective potential of the dynamical mean-field theory, defined 
by 



e -v DMFT [™+] = . qC)C uexp< - > n int [c ia ,c. 



/ d[c, c + ] exp <^ - H h 

ka ) 



(10) 



Gok = Gok — C{iVn) is the bare Green function of the non-local degrees of freedom. 
Similarly to V[r],r] + ], the functional Vdmft[?7, V + ] generates connected vertices (which 
are in general one-particle reducible), amputated by the bath Green function £(ii/ n ). 
To simplify Eq. (jHJ), we perform a shift Cka — > Ck a — Vka, such that 

, / rr~~+l f^--- r 'ta)Gok(^-Vk^)-VBMFT[c,C+] 

V[?7,?7 + ] = -In / d[c, c + ]e fc ' CT . (11) 



To arrive at the standard dual fermion approach [16j [T71EIE1 [19] we consider an expansion 
of Vdmft [c,c + ] in fields 

V DM ft[c,c+] = _ alZZi(iuy ^ + ^E^l 



k.cr 



V D mft[c, c + ] = ^r loc o (4 iCT C fe2 J(c+ 3(T ,c fe4CT (12) 

+ l^loc ° (Cfe 1 /fe 2 J(c£ (T ,C fc4CT 0(cJ (T ,,C fc6CT ») + ... 
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where ri oc and T^ c are the connected 4- and 6- point local vertices, amputated with the 
bare Green functions (, e.g. 

rgii^M) = (i + <wr x II 4 , C~V<) (13) 

x ^loc.aa-'C^i-^s) - G , i oc (iz/i)Gi oc (iz/ 2 )((5 I/li/ 3 - <W<W 3 ) 

and o stands for summation over momenta- frequency- and spin indices fulfilling the 
conservation laws, G[^ c is the two-particle local Green function, which can be obtained 



via the solution of the impurity problem|10[ [22]. We therefore obtain 

v[v,v + ] = -E<A*V* (14) 



k.c 



In / d[c, c + ]e k ' a 



E [%A G 'k) lc k,<y~K,a G Ok 1k, a -VZ,a G 0k C fc,aj-V D MFT[C, 



■ 



where 



G' k = gJ—^P^^ = [1 - C(^n)S loc (^ n )] 2 G k , (15) 

and Gk = Gk — G\ QC (iuj n ). Rescaling the fields of integration to exclude extra factor 
(1 — CHioc) 2 and introducing the 'dual' source field 

rjka = Vka/ [1 ~ Si oc (w n )G fc] , (16) 

we obtain the effective interaction of the lattice theory in the form 

V[ V ,r] + ) = V[%rt\ ~ - 1 G, 1 ) ^ (17) 



1 - Zj| n „ l^n Uflt 



fc,<7 

where 



^loc 



E G fc 1 (c+ -»?+ ) (c fc CT -?? fe a ) - Vdmft [c + ,c] 

V[r7,7^] = -ln / D[c,c + ]e fe - V ' ' A ; (19) 



is the effective interaction of dual fermions. According to the Eq. ( TT2|) . the expansion 
of Vdmft [c + ,c] in fields reads 

Vdmft[c + ,c] = ^r ioc o (c+ ua c k3 ,a)(c+ 2 a ,Ck 4 y) (20) 

+ l T iol ° (cjfc 1 /fe 2 J(c^ 3CT ,c fc4(J /)(c+ <T ,,c fe6(J ») + ... 

where r ioc = f ioc nt=i(! ~ C s ioc)i^ and = f [IiUC 1 _ C s ioc)^ are the connected 
4- and 6-point vertices, amputated with the local Green functions G\ oc . For the four- 
point vertex ri oc the requirement of connectivity and amputation with the full local 
Green functions implies one-particle irreducibility. However, the higher-order vertices, 
e.g. r{^2 remain one-particle reducible. 



a) b) 

• *-•> • • • 

Figure 1. (Color online) (a) The one-particle reducible contributions to the lattice 
two-point Green function, coming from one-particle reducible six-point local vertex 
(shaded). The circles correspond to the four-point local vertices, solid line - the 
propagator Gfc, dashed line corresponds to the local propagator G/ oc . (b) The 
corresponding contribution to the non-local self-energy, which is of the second order in 
the local 4-point vertices. 



In the described approach the relation between the self-energies £d(&) of the 
dual fermions c& and £& of the lattice fermions is easily obtained by equating the 
corresponding two-point vertices in Eq. ( ITT)) : 

Sfc _ £i oc (zz/ n ) 
1 — S fc Gofe 1 — Yii oc (ii> n )Gok 

+ I (21) 

l-S d (fc)[G fc -G loc (^ n )][l-S loc (z^)G 0fc ] 2 ' 1 ' 

which implies 

S fc = + E loc (iv n ). (22) 

The result (122)) was derived in Ref. [T6] . 



3. The effect of the six-point vertex 

3.1. Relation between the dual and lattice self-energy 

The relation (1221 does not change its form in the approximation, when one neglects six- 
point (and higher order) local vertices in Eq. ( fT2l) . This however does not necessarily 
mean, that it remains correct in this case. Instead, as becomes clear from the following 
discussion, the relation (122)) implicitly assumes that the one-particle reducible diagrams 
for six- and higher order vertices are taken into account. 

Let us consider the diagrams for the self-energy, which include six-point one particle 
reducible vertex, such as the diagram shown in Fig. la. These diagrams, being 1PI in 
terms of the vertices , produce nevertheless one-particle reducible contributions to 
the self-energy, which should be excluded. The denominator in the relation ( |22l) aims to 
remove the corresponding one-particle reducible diagrams for self-energy. Formulated 
differently, the quantity £' in Eq. (122)) must contain one-particle reducible diagrams, 
which are cancelled by the denominator in the first term in Eq. ( 1221) . 

To prove this statement for the diagrams, similar to that of Fig. la, containing 
repeating lowest (second-order) diagram of Fig. lb, it is sufficient to consider the 
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tree diagram contribution to the six-point vertex of the form = ^2 Fi oc Gi oc ri oc , 
the sum is taken over different combinations of 4-momenta, assigned to 4-point 
vertices. Decoupling the resulting six-particle interaction with the fermionic Hubbard- 
Stratonovich transformation by introducing auxiliary fermionic field <pk, we obtain 



e -V[rf,V+,C,C + ] 



= J d[c, C + ] exp I [ G k 1 ( C k,a - Vk,a) i C k,a ~ Vk,a) 

+^(^n)(0+ CT - C+Mk* - Ck,a)] (23) 

+ Tioc o (cJ j<T c &3iff )(c'+ )(T ,0 fe4i(r /) 

+ ri oc o (c+ iCT c fc3i(T )(0+ i(T ,c fc4i(T /) (24) 

where we have introduced source fields ( for fermions <fi. The effective interaction 
(|23|) can be put in more compact form by introducing spinors $fc >cr = (c^, 4>k,a) and 
n fcCT = (r) k ,a, (k,<r), such that 

e -vp,n+] = / d[$,$+] e fc .- (25) 



where the corresponding matrix bare Green function reads 

Gk = ( dh ° x 1 (26) 
y G\ oc [iv n ) J 

It is of crucial importance that the matrix Green function, Eq. (1261) . contains both, the 
non-local G k and the local G\ oc (iu n ) components, which are mixed through the non-local 
dual self-energy, as considered below. 

The resulting two-point vertices can be also considered as matrices in the space 
(fj, Q. The relation between the lattice and dual two-point vertices then has the form, 
similar to the Eq. ( fTTl) . 

t Sloc(Wn) n \ 

i-s loc (w n )G 0fe u + yW (27) 
J k 

V k {2) =Z d (k)ll-G k Z d (k)]-\ 

where E d (fc) is the self-energy matrix of c k ,a and 0fc jCr fields, having both, diagonal and 
off-diagonal contributions. For the second-order diagram of Fig. lb the self-energy is 
equal for both fermion species: 

S d (*) = Ej^fcW j M. (28) 

where E d PI (fc) is the value of the IPI diagram of Fig. lb. From this we obtain the 
relation between the local and lattice self-energies 

r (2)c,c _ S loc (^ n ) £ d PI (fc) 1 



1 - 1 - E Ioc (i^)G , 1 - G k Zf\k) [1 _ E loc ^ n )G ofc f 
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a) b) 




Figure 2. (Color online) The non-local contributions to the lattice self-energy, which 
are of the third order in four-point local vertices, (a) The contribution, containing 
the local Green function, coming from the one-particle reducible six-point vertex, 
(b) Similar contribution of the dual fermion approach, neglecting six-point vertices. 
Notations are the same, as in Fig. 1 



= ; * , (29) 

1 — Z^Crofc 

which yields 

S fe = S loc ( Wn ) + ^ PI (^)- (30) 

The result ([30]) is essentially different from Eq. (122]) and implies that the one-particle 
reducible contributions in the self-energy, occurring due to the one-particle reducible 
contributions to the six-point vertex, are indeed cancelled by the denominator in Eq. 
(122]) . Equtions ([221 and ( [30]) also imply the analogue of the Dyson equation for the dual 
fermion self-energy 

SdW = l-GuJjEFikY (31) 

Having the IPI self-energy S^ PI (/c) the result (I3T)]) should be used to obtain the lattice 
self-energy instead of the equation (!22j) . suggested by Refs. [HI [T7J HH [19]. 

For the considered theory with only six-point one-particle reducible contributions 
included, the newly derived relation (13"U]) between the dual and lattice self-energy is 
fulfilled if (and only if) the self-energy is equal for both fermion species c and </>, as 
it happens for the diagram of Fig. lb. The abovementioned assumption does not 
necessarily hold in higher orders of dual perturbation theory. However, the inclusion 
of one-particle reducible contributions of higher order (eight and more point vertices) 
makes the relations (I3(jp and (1311) fulfilled in more general situations[23j. 

The one-particle reducible contributions to the six-point vertex can also produce IPI 
self-energy diagrams, containing local Green functions, such as shown in Fig. 2a. These 
diagrams, being formally of the same order in the local 4-point vertices, as considered 
by the dual fermion approach (Fig. 2b), are not taken into account when the six- 
point and higher vertices are not taken into account. At the same time, the diagrams, 
similar to that shown in Fig. 2a, can produce even larger contribution to the self-energy, 
than the diagrams of the dual fermion approach, neglecting 6-point and higher vertices, 
since the sum of the local Green function over momentum does not vanish. Therefore, 
accounting one-particle reducible parts of six- and higher-order local vertices appears 
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to be crucially important for both, the diagrammatic consistency of the dual fermion 
approach and keeping all the diagrams of the same order in the four-point local vertices. 

We also note that the covariation splitting method, used in the dual fermion 
approach (19~1), is similar to that applied in functional renormalization-group approach 
(see, e.g., Refs. [2U El]), except that the latter considers integration of degrees of 
freedom in many infinitesimally small steps, while the former - only in two steps. 
Similarly to the discussion above, in the Polchinski formulation of the functional 
renormalization-group approach [24"l [21] one-particle reducible contributions to six-point 
vertices were argued to be important for proper calculation of the four-point vertices 
already in one-loop approximation [25]. The same contributions can be also shown to 
be important for evaluation of multi-loop contributions to the self-energy. In the next 
subsection we address another aspect, where the six- and higher-order vertices appear 
to be important in the dual fermion approach. 

3.2. Self-energy feedback 

The dual fermion approach accounts for the self-energy feedback by dressing the Green 
function of the non-local degrees of freedom: 

G k = G k /[l-G k E d (k)}. (32) 

However, the function G k , which represents a difference of two propagators, does not 
correspond to a physically observable quantity, and it is informative to trace, how 
dressing it one can finally obtain the physical propagator 

G k = G k /[1 — G k [E k — Sioc)] (33) 
= G k /[l-G k Z d p \k)}, 

which is constructed by dressing the Green function G k , see Eq. ([6]), containing only the 
local self-energy, by the remaining self-energy difference — Si oc . We have observed in 
Sect. 3.1, that in the lowest orders of perturbation theory, propagators G k , appearing 
in the diagram technique for Eq. ffl9|) . are added by G\ oc , coming either from either 
adding local quantities to their non-local counterparts (such as in Eq. ( |30|) ). or from 
the contributions, containing one-particle reducible six-point and higher order vertices, 
such as the diagram of Fig. 2a. Adding G\ oc to G k is however still not sufficient to 
reproduce ([33]) for Z\ pl (k) ^ 0. 

Again, we argue, that the six- and higher order local vertices are crucially important 
to obtain ( 1331) . To see this, let us insert G k = G k + Gi oc into Eq.( l33|) . use Eqs. ( I3TI) 
and (|32|) to represent the result in terms of G k and H d PI (k), and expand the result in a 
series of T} pl {k) : 

CO 

G k = Yptoc&n) + nG k ][G loc (iv n )Zf '(k)]^ 1 (34) 

n=l 

The term with n = 1 represents the combination G\ oc + G k , discussed above, while 
the terms with n > 2 in this series expansion can be ascribed to the respective diagrams 




Figure 3. (Color online) Dressing of the non-local Green function G^, shown by bold 
line and appearing as a part of diagram (a) by the two contributions of the lowest-order 
in E^ PI (fc) (b). Notations are the same, as in Fig. 1 

(see Fig. 3 for n = 2), where the one-particle irreducible local vertices are connected by 
local propagators, forming one-particle reducible six-point and higher vertices. The dual 
fermion approach, which does not account for the six-point and higher order vertices, 
neglects therefore a difference between and + G\ oc . 

4. Conclusion 

In the present paper we have considered effect of one-particle reducible six- and 
higher-point vertices in the dual fermion approach. We have argued that the one- 
particle reducible contributions to these vertices are important to make the dual 
fermion approach diagrammatically consistent. Neglecting the six-point and higher 
order vertices does not allow to obtain correct relation between the dual and lattice 
self-energies, as well as treat correctly the feedback of the dual self-energy on the dual 
Green functions. Apart from that, the one-particle reducible six-point and higher order 
vertices lead to the self-energy corrections, which contain both, local and non-local 
Green functions. 

Further numerical investigations of the (un) importance of the described 
contributions of six-point and higher-order vertices to be performed (see, e.g., Ref. 
[23J). This also calls for developing a new method of treatment of non-local degrees of 
freedom, which avoids the described problems of the dual fermion approach and operates 
with the one-particle irreducible quantities. 
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